Helical vortex phase in the non-centrosymmetric CePtaSi 
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We consider the role of magnetic fields on the broken inversion superconductor CePtsSi. We show 
that upper critical field for a field along the c-axis exhibits a much weaker paramagnetic effect than 
for a field applied perpendicular to the c-axis. The in-plane paramagnetic effect is strongly reduced 
by the appearance of helical structure in the order parameter. We find that to get good agreement 
between theory and recent experimental measurements of Hc2, this helical structure is required. We 
propose a Josephson junction experiment that can be used to detect this helical order. In particular, 
we predict that Josephson current will exhibit a magnetic interference pattern for a magnetic field 
applied perpendicular to the junction normal. We also discuss unusual magnetic effects associated 
with the helical order. 
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PACS numbers: 

The recently discovered heavy fermion superconductor 
CePtsSii has triggered many experimental and theoreti- 
cal studiesSiSi^iSiSiLSiS. There are two features which have 
caused this attention: the absence of inversion symmetry; 
and the comparatively high upper critical magnetic field 
{Hc2)- Broken inversion symmetry (parity) has a pro- 
nounced effect on the quasiparticle states through the 
splitting of the two spin degenerate bands. This influ- 
ences the superconducting phase, which usually relies on 
the formation of pairs of electrons in degenerate quasi- 
particle states with opposite momentum. The availabil- 
ity of such quasiparticle states is usually guaranteed by 
time reversal and inversion symmetries (parity It 
is relatively easy to remove time reversal symmetry, e.g. 
by a magnetic field, and the physical consequences of 
this have been well studied. However, parity is not so 
straightforwardly manipulated by external fields. Super- 
conductivity in materials without inversion center there- 
fore provides a unique opportunity in this respect. 

The large Hc2 ~ 4r in CePtaSii*^ imphes that the Zee- 
man splitting must be non-negligible below Tc — 0.75K 
(the estimated paramagnetic limit is at Hp w 1.2 T). In 
a magnetic field, this superconductor has to form Cooper 
pairs under rather odd circumstances. In particular, it 
is no longer guaranteed that a state with momentum k 
at the Fermi surface has a degenerate partner at —k. 
The state k would rather pair with a degenerate state 
—k + q and in this way generate an inhomogeneous su- 
perconducting phase. We argue below that recent Hc2 
measurements^ suggest that this is the case in CePtaSi. 
These measurements show that, while the upper critical 
field is basically isotropic close to Tc, a small anisotropy 
appears at lower temperature^ {H^2/H^2 = 1-18 at 
T — Q). The apparent absence of a paramagnetic limit in 
CePtaSi can be explained by lack of inversion symmetry 
even if the pairing has s-wave symmetrjiSiiSii^. However, 
these works indicate that suppression of paramagnetism 
is very anisotropic and the application of this theory to 
CePtaSi would indicate no paramagnetic suppression for 
the field along the c-axis, but a suppression for the field 
in-plane (ffp^ w 1.7T). The relative lack of anisotropy 



in the experimental data is surprising in this context. 

In this letter we examine the mixed phase close to 
the upper critical field. Using the crystal symmetry of 
CePtaSi, we show that the high-field superconducting 
phase has pronounced differences for field-directions par- 
allel and perpendicular to the c-axis. For the field parallel 
to the c-axis, the paramagnetic limiting is basically ab- 
sent and the vortex phase is quite conventional. While 
for the perpendicular field direction, the field can induce 
a phase which gives rise to an additional phase factor in 
the many body wavefunction e*'^' -^ with q perpendicular 
to the applied field: a helical vortex phase. We also pro- 
pose a Josephson junction experiment that can be used 
to detect this helical phase factor and discuss a transverse 
magnetization related to the helical phase. 

We use the single particle Hamiltonian 



H 



^'^ks' (1) 



where (cfe^) creates (annihilates) an electron with mo- 
mentum k and spin s. The band energy Cfe = ^fc ^ 
measured relative to the chemical potential /i, agj^ ■ cr 
introduces the antisymmetric spin-orbit coupling with a 
as a coupling constant (we set = 1 where (..) is an 
average over the Fermi surface), and ■ cr gives the 

Zeeman coupling. The crucial term in Eq. ^is agj^ ■ cr, 
which is only permitted when inversion symmetry is bro- 
ken (g^ satisfies gj^ — —g_]f due to time reversal symme- 
try). This term destroys the usual two- fold spin degen- 
eracy of the bands by splitting the band into two spin- 
dependent parts with energies Ej^ ± = ^fc '^ISfcl- 
spinors are determined by the orientation of the corre- 
sponding gj^ . The general pairing interaction is 



T~(-pair — 2 ^ ^ ^(^7 ^ ) 

k,k,q,si 



^ 4:+q/2.si'^-k+q/2.s2^-k'+q/2,s2^k'+q/2,si' 

(2) 

expressed in the usual spin basis. We will work in the 
large a limit so that the pairing problem becomes a real 
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two-band problem in the diagonal spinor (±) basis. To 
find the pairing interaction in the ± basis, we diagonalize 
the single particle Hamiltonian after which the two-band 
pairing interaction, for H = 0, is written in spinor form 
as 

y=-mfc')(:.,4: u^x) (3) 

where 0± = 0^ ± (f>j^i and A± — (1 ± • g^/) where we 
have taken g/, = \gf, \ ( sin 0/, cos 0^ , sin Oj^ sin (j)f^ , cos 0^) . 
Note that even for a spatially isotropic interaction, the 



two-band solution has both a spin-triplet and a spin- 
singlet gap function when a (this is a consequence 
of the broken parity symmetrjiii) . We will consider the 
limit a » ^bH and keep only terms up to order iisH/a 
(a good approximation for CePtaSi). We restrict our- 
selves to choices of ^(fc, fc') that corresponds to spin- 
singlet pairing in the a = limit. This restriction allows 
us to use Eq. 13 even if ff 7^ 0, which considerably sim- 
plifies the notation. 

With the two band pairing interaction of Eq. O the 
linearized gap equation becomes 



A„(fe, q)=-Tj2Yl + 9/2. i^n)G${-k' + q/2, -zu;„)A^(fc', q) (4) 

n,k' 

I 



where Gj-{k,iujn) = [iuJn + 5fc ^ c^lfl'fel] ^- 
ting A+(fc,g) = e'^kA+{k,q) and A_(fc,g) = 
— e~*'^fc A_(fe, q) results in the simplified two-band equa- 
tion with the interaction 

y-y(fe,fe')(fTo-tT,)/2. (5) 

The factors g in Eq. do not appear because of Pauli 
exclusion and the assumed singlet form of V{k, k') in the 



a = limit. We denote the density of states on the 
two Fermi surfaces by iV+ = Ncos^{5/2) and N_ = 
Ns\-a^{5/2). We can write Ka{k,q) = ?Ar,Q(g)/r(fe), 
where Eq. [S] implies that ■!/'r,+ ('Z) + V'r.-(9) — 0. With 
Vr.+ lg) = — V'r,-('Z) = ^{q) and the proper Fourier 
transform of the gap equation keeping gauge invariance, 
we find the following equation determining the upper crit- 
ical field 



*(H)i^= / . '^^ (\frik)\^[cos{pH ■gu) + i cosS si 
t Jq smh(mj \ 



where t = T/T,, H = Hvl/{27r<i>oT^), jl = /is2r,$o/t'F, 
^ (I'D = {vl + v± = {vi ± iv2)/\/2, vi^2 are 
components of the Fermi velocity perpendicular to the 
magnetic field, n± = {Di^iD2)lH / V^, In = ^o/(27ri7), 
and Di — —idi — 2eAi/c. The upper critical field Hc2 is 
found by expanding ^[R) — J2n^ri4>n{R) {4>n{R) are 
the usual Landau levels). 

In the following we take a spherical Fermi surface. For 
CePtaSi this will not be the case, but the overall geome- 
try of the Fermi surface does not qualitatively change our 
results. We also take gj^ = ^/3/2{-ky, k^, 0) as the lowest 
order term in k allowed by symmetry and consider the 
case V{k, k') = Vo for isotropic s-wave pairing. Our re- 
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suits will hold for any pairing symmetry as the Ginzburg 
Landau (GL) theory discussed later will demonstrate. 

For the field along the c-axis, g ■ H = so Hc2 is 
independent of the Zeeman field and there is no param- 
agnetic effect (note that in principle there can be a para- 
magnetic effect since there are allowed by symmetry 
that contain a component^, however such terms are 
expected to be small) . The solution of the upper critical 
field problem is identical to that carried out by Helfand 
and Werthamer^** and is plotted in Fig. 1. However, for 
fields perpendicular to the c-axis unusual properties oc- 
cur, which can be best illustrated by a GL theory with 
free energy density 



/ = -a|*|2 + q^-l^ + — + — en • B X [*(£»*)* + **(£>*)] + — (6) 

2 Im ZlTlc OTT 

I 



where a — ao(T° — T), D — [D^, Dy), n is the unit vector oriented along the c-axis, and B = V x A. Eq.|Slapplies 
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FIG. 1: Upper critical fields for CePtaSi with fields along 
the c-axis and in the plane. The actual in-plane Hc2 will lie 
between the two extremes shown. The efi'ect of the helical 
order on Hc2 can be quite pronounced. These calculations 
are for aM =3. 

to all possible pairing symmetries with a single complex 
order parameter, as discussed also by Samokhinl. The 
lack of inversion symmetry allows for the existence of 
the term proportional to e (for a discussion of other re- 
lated terms see Ref . Iisj) . This term induces a spatially 
modulated solution in a uniform magnetic field. The GL 
equation for the order parameter ^{R) = e^1'-^'^{R) in 
a field is identical to the zero-e GL equation with 'l'(fi), 
where 

q = -2me n x B (7) 

(a microscopic expressions for q is given below). Conse- 
quently, the upper critical field solution in the GL limit is 
"^{R) = 4>Q{R)e^^'-^. We call this phase the helical vor- 
tex phase. The helical order coincides with an increase 
in the upper critical field {B = H): 

T,(H)=T,- + \ ^. (8) 

The expression for the supercurrent density is 

J = cV X [B - 47rm] /47r = Jo + 4ee(n x (9) 

where J q is the usual supercurrent density in the GL 
limit and m — x Jq is the magnetization due to the 
e term Eq. |^ The usual boundary condition on the order 
parameter is given by J = through the boundary. The 
appearance of m in Eq. |51is highly unusual in GL theory 
and has some consequences that are discussed later. Note 
that if l^*! and B are spatially uniform then J = for 
q given by Eq. [7| the helical phase carries no current. 
The possibility of this helical order has been raised in 
the context of thin film or interface superconductivity 
where the vector potential can be neglectediSiiLiS. As 
discussed below, we find here that it can play a very 
important role in the vortex phase. 

The increase in Hc2 due to the appearance of the helical 
plays an important role in the microscopic theory. Since 



q- 0, we can expand (^o(-R)e'^ = K{q)4'n{R) 
where 6„ = {iqlny^e'^'^^"'' 1^ l\/2^n\. Our numerical mi- 
croscopic solution has this form and near T^, we find, 

1 - 2,bHcosS J^;^TJJ^^ , (10) 

for /r(fc) = 1 this gives qln — Q.'iTiaM cos 5 \/WJh^ 
where um = ^J^H^il Hp is the Maki parameter, li^^ 
is the upper critical field for /is = (this coincides 
with the i/c2 for the field along the c-axis), and Rp — 
A/(/XB\/2). Fo r /r(fc) = kl - Eq.UHl gives qIh = 

0.418Q;Af cos Sy^ H/H^2- The enhancement of Hc2 due to 
the helical order can be substantial as Fig. 1 shows. Our 
numerical results show qln can be larger than one, which 
implies that the helical wavelength becomes less than the 
spacing between vortices. Fig. 1 is for aM — 3 which is 
slightly smaller than um — 3.8 that follows from the 
measurements of Ref. 0- The helical order changes with 
varying cos 5. For cos (5 = the density of states are 
the same on both Fermi surfaces and no helical order 
appears (labelled 'in-plane (no helical)' in Fig.l); while 
for cos (5 = 1, qln is maximum (this corresponds to the 
curve labelled 'in-plane (helical)' in Fig. 1). For all other 
possible values of cos 6, the Hc2 curve lies between these 
two extremes. The limit cos (5 = 1 is unlikely since this 
implies that the density of states of one of the two bands 
vanishes. However, if the elements of the pairing interac- 
tion Va^f3 in Eq. are different in magnitude from each 
other, then a large qln can still arise for cos (5 = 0. 

Comparing our result with the i7c2-ineasurement by 
Yasuda et al.^ for CePtaSi, we find that only taking 
the effect of the spin-orbit coupling into account would 
not account for the relatively large value of the in-plane 
Hc2 ~ 2.8T) at T 0. Here, paramagnetic limiting 
should reduce the value to below 2 T. We can, however, 
explain the increased Hc2 by including the helicity. 

We have also examined the behavior of the Abrikosov 
parameter Pa = (l^(-R)l^)/(l*(-R)P)^ in connection 
with the possible vortex lattice structures. We find a 
possible structural transition from a stretched hexagonal 
lattice at high temperatures to a stretched square lat- 
tice at low temperatures. The origin of this transition 
is related to the two-dimensional inhomogeneous state 
discussed in Refs. IJs, 19 20. Note that the helical order 
discussed above is distinct from that discussed in these 
works. The physics discussed in Refs. 118^19. does not play 
a significant role in CePtaSi because the value of um is 
too small. However, if Q!a/ >> 1 then the vortex physics 
becomes very exotic. 

We have not discussed the direct experimental verifi- 
cation of the helical phase. Since helicity of the order pa- 
rameter is related to its phase, an interference experiment 
based on the Josephson effect would provide the most 
reliable test. Here we propose to consider a Josephson 
junction between two thin film superconductors (Fig.[5J), 
one (1) with and the other (2) without inversion symme- 
try (for CePtaSi the c-axis is perpendicular to the film). 



4 



Inversion 
symmetric 



H 



Broken 
inversion 




FIG. 2: Josephson junction geometry for the observation of a 
helical phase. 



For a magnetic field applied in the plane of the film per- 
pendicular to the junction and with the superconductor 
(2) is oriented so that the helicity q is perpendicular to 
the field ; we find this gives rise to an interference effect 
analogous to the standard Fraunhofer pattern. For this 
experiment, the film must be sufficiently thin that the 
magnetic field and the magnitude of the order parameter 
are spatially uniform. 

To illustrate this, consider the following free energy of 
the junction 



Hj = -t / dx[^i{R)^l{R) + c.c] 



fll) 



where the integral is along the junction. The resulting 
Josephson current is 



Ij^Im t dx^i{R)^2{R) 



(12) 



Setting the junction length equal to 2L, and integrating 
yields a maximum Josephson current of 



/; = 2i|*?||vl>0| 



|sin(g-L)| 



(13) 



This, combined with the result of the microscopic the- 
ory near (for an isotropic interaction) that qL — 
2.4:COsSaMH^oL/^o (where = 0-lSvp/Tc) demon- 
strates that the Josephson current will display an inter- 
ference pattern for a field perpendicular to the junction. 
Note that in the usual case the Fraunhofer pattern would 
be observed for a magnetic field perpendicular to the thin 
film. Furthermore, for a sufficiently large in-plane field, 



so that \qL\ >> 1 (which implies /j « 0), the arguments 
of Ref. |21| imply that the application of an additional 
magnetic field along the surface normal will lead to an 
asymmetric Fraunhofer pattern. 

A less direct probe of the helical order is to look 
for the related transverse magnetization that appears in 
Eq. O There are two situations for which this can be 
observed. The first is in a thin film with a supercur- 
rent flowing along x and a surface normal along y. In 
this case a magnetization will exist along y (normal to 
the film). This situation is a generalization of an ex- 
periment originally proposed by EdelsteinpSS. We esti- 
mate \m\ = ^ cosSusfiBVs/vF ~ 0.02 Gauss assuming 



Vs/vf = 2 X 10 cos (5 — 1/3, and 



1 X 10^8 m-3. 



The second is in the vortex lattice phase for a field ap- 
plied along the c-axis. In this calculation valid 
near Hc2 gives m — x Jq, where Jo is the usual 
supercurrent for the Abrikosov vortex lattice. Near the 
vortex core, m is directed radially outward, perpendicu- 
lar to the applied field. 

We have considered the role of magnetic fields on the 
non-centrosymmetric superconductor CePtsSi. Using a 
two-band theory with a Rashba spin-orbit interaction, 
we have shown that the upper critical field for the field 
along the c-axis behaves as if it would in a conventional 
superconductor independent of the paramagnetic (Zee- 
man) field. We have further shown that while there is a 
paramagnetic limiting effect for magnetic fields applied in 
the basal plane, this effect can be strongly reduced by the 
appearance of a helical order. Our theory agrees with the 
experimental measurements of i?c2, despite a relatively 
strong Zeeman field, provided this helical order exists. 
Finally, we have proposed a Josephson junction experi- 
ment that can unambiguously identify the helical order 
and discussed the appearance of a transverse magnetiza- 
tion related to the helical order. 
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